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Abstract. We introduce a new method for computing heat invariants 
a n (x) of a 2-dimensional Riemannian manifold based on a commutator 
formula derived by S.Agmon and Y.Kannai. Two explicit expressions 
for a„(x) are presented. The first one depends on the choice of a certain 
coordinate system; the second involves only invariant terms but has some 
restrictions on its validity, though in a "generic" case it is well-defined. 

1. Introduction and main results 

Let M be a (i-dimensional compact Riemannian manifold without bound- 
ary with a metric {gij}- Let A be the Laplace-Beltrami operator (or simply 
the Laplacian) on M. In the coordinate chart {u{\ it is written in the form 

where g =det{gij}, and {g lJ } denotes the inverse of the matrix {gij}. 

The Laplacian is a self-adjoint elliptic operator acting on smooth functions 
/ : M — ► R. Let {Aj, be the spectral decomposition of A (considered 

as an operator in L 2 (M)) into a complete orthonormal basis of eigenfunc- 
tions ifii and eigenvalues Aj, < Ai < A2.... Consider the heat operator e~ tA 
which is the solution of the heat equation (d/dt + A)f = 0. For t positive, it 
is an infinitely smoothing operator from L 2 (M) into itself and has a kernel 
function 

K(t, x 1 ,x 2 ) = e~ Xlt (fi(xi) (g> <pi(x 2 ), 
analytic in t and C°° in x\ and x 2 , such that 

(e"* A /)(xi) = J K{t,x l ,x 2 )f(x 2 )dx 2 . 

M 

The following asymptotic expansion holds for the trace of the heat kernel: 



00 



d 



K(t, x,x) ~ 2^ a n (x)t n 2 

n=0 

as t — > + ([Gl]). The coefficients a n {x) are called heat invariants of the 
manifold M; they depend only on the germs of the metric. Moreover, it is 
known that the heat invariants are homogeneous polynomials of degree 2n 
in the derivatives of {g lJ } at the point x ([G2]). 
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Integrating a n {x) over the manifold one gets the coefficients of the expan- 
sion of the trace of the heat operator: 



where p(x)dx is the volume form on the manifold M. 

Heat invariants are metric invariants of the given manifold M and contain 
a lot of information about its geometry and topology (see [Be], [G2], [G3], 
etc.). For example, ao determines the volume of M, a\ its scalar curvature 
(and hence the Euler characteristic in dimension 2). Higher coefficients are 
of a considerable interest to physicists since they are connected with many 
notions of quantum gravity (see [F]). 

The structure of the heat invariants becomes more and more complicated 
as n grows; compare the formula for 03 ([Sa], 1971) with the formulas for 
a 4 ([ABC], 1989; [Av], 1990) and a 5 ([vdV], 1997). For n > 5 no explicit 
expression for the heat invariants is known, even though there are many 
algorithms to derive them — the computational complexity is so high that 
even with the help of modern computers it is impossible to get the result. 

In the present paper we introduce a new method for computing heat 
invariants. We call it the commutator method since it is substantially based 
on a commutator formula by S. Agmon and Y. Kannai ([AK]) which deals 
with asymptotic expansions of the resolvent kernels of elliptic operators (let 
us note that the asymptotics of the relevant kernels associated with a self- 
adjoint realization of the elliptic operator depend only on the local behavior 
of its coefficients). 

In section 2 we obtain the following concise reformulation of the Agmon- 
Kannai formula: 

Theorem 1.2. Let H be a a self-adjoint elliptic differential operator of 
order s on a manifold M of dimension d < s, and Hq be the operator 
obtained by freezing the coefficients of the principal part H' of the operator 
H at some point x G M: Hq = H' x . Denote by R\(xi,X2) the kernel of 
the resolvent R\ = (H — A)" 1 , and by F\(x\,X2) — the kernel of F\ = 
(Ho — A) -1 . Then for A — > 00 there is a following asymptotic representation 
on the diagonal: 




00 



R x (x,x) ~ p(x)' 1 X X m F™ +1 (x,x), 



(1.3) 



m=0 



and the operators X m are defined inductively by: 



Ao — I-, X m — X m _iHo — HX m -\. 



The commutator method may be applied to various problems of "heat 
kernel"-type (one may consider other operators instead of Laplacian, Lapla- 
cian on forms, etc. — see also [Kan] for the Agmon-Kannai formula in case 
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of matrix operators). In this paper we are focusing on the simplest non- 
trivial question — what are the heat invariants of a 2-dimensional Riemann- 
ian manifold? 

An advantage of the commutator method is that it yields very explicit 
formulas for all 2-dimensional heat invariants (cf. [Xu] where heat invariants 
of manifolds of arbitrary dimension are considered) . They are given by the 
following 

Theorem 1.4. Let (u,v) be local coordinates in a neighborhood of the point 
x = (0, 0) on the Riemannian 2-manifold M, in which the metric is confor- 
mal: 

ds 2 = p(u, v)(du 2 + dv 2 ). 
Then the heat invariants a n (x) are given by 

An m k—n 

«n(*)= E E J2 C nksmP(u,v) k - n A k (u 2k - 2n - 2s V 2s )\ u=0 , v=0 , 

m=n+l k=n+l s=0 (1-5) 

where C n k S m are constants given by (3.12). 

The proof of this theorem is presented in section 3. We also perform a 
test of the formula (1.5) which is described in the Appendix. 

Let us note that local coordinates in which the metric is conformal always 
exist (see [DNF]). We also observe that the derivatives of p(u,v) appear in 
(1.5) through the formula (1.1) for the coefficients of the Laplacian. 

The expression (1.5) depends on a coordinate representation but it is 
possible to make it invariant rewriting it in special curvature coordinates 
which are introduced in section 4. The resulting formula is no longer valid 
in general but is true "generically" . 
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under the supervision of Professor Y. Kannai. I am very grateful to Profes- 
sor Kannai for his constant help and encouragement. I would also like to 
thank Professor L. Polterovich and Professor M. Solomyak for many fruitful 
discussions, and my colleagues J. Greenstein and A. Grigoriev for helpful 
advice. 

2. Multiple commutators 

We recall the notion of a multiple commutator of two operators, intro- 
duced in ([AK]), which plays a fundamental role in the sequel. 

Definition 2.1. Let J be a finite vector with nonnegative integer compo- 
nents and let A, B be linear operators on some linear space M. The multiple 
commutator [B, A; J] is inductively defined for all such vectors J in the fol- 
lowing way: 

1. [B,A;0]=B 
2. [B,A;j] = [[B,A;j-l],A] 
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3. [B,A;JUj] = [B[B,A,J],A;j], 

where J U j denotes the vector obtained by adding the component j to the 
vector J to the right. 

We introduce a filtration on the space of all vectors with non-negative 
integer entries. Denote by V m = {J = (ji,-,j r ) '■ \J\ + r = to}, where 
\J\ =3i + ■■+ Or- Let 

X m = J2 [B,A;J], rn>\. (2.2) 

JeVm 

Theorem 2.3. The operators X m satisfy the following recurrent relation, 
X l = B; X m = BX m ^ + [X m -i,A], (2.4) 
and are given explicitly by 

m \{A-B) k A m - k . (2.5) 

Proof. Let us proceed by induction. Indeed, X\ = B since V\ consists of 
a single vector (0). In order to pass from X m -\ to X m we do the following. 
Let denote the set of all vectors in V m whose last entry is 0. Consider a 
mapping p : V m — ► V m -\ defined as follows: for J = (ji, ..,j r -i,0) G we 
have p(J) = (ji, ..,j r -i) and for J = (j 1 ,..,j r ) G V m \ V,® n we have p(J) = 
{hi—>3r — !)• It is clear that the mapping p is well-defined. Moreover, 
it is a surjection and every element of V m -i has exactly two preimages - 
one in V^ n and one in its complement. Applying the definition of a multiple 
commutator we prove the inductive step — BX m _i is the contribution of 
vectors from and [X m -i, A] is the contribution of vectors belonging to 
V rn \ V^. This completes the proof of the first part of the theorem. For 
obtaining the closed formula, let us rewrite the recurrent relations in the 
following way: X\ = A - (A - B) = B; X m = X m ^A -(A- B)X m _ 1 . 
Proceeding by induction and recalling that (™) + ( fc ™ x ) = ( m ^J _1 ) we get 
(2.5) ' " □ 

Let us present the original Agmon-Kannai formula (see [AK]): 

Theorem 2.6. In the notations and conditions of the Theorem 1.2 the fol- 
lowing asymptotic representation on the diagonal holds for the kernel of the 
resolvent R\ : 

oo 

R x (x, x) ~ p(x)-\F x (x, X ) + J2([H - H, H ; J]F^+ r+ \x, x))), 

J (2.7) 

where the sum is taken over all vectors J of length > 1 with nonnegative 
integer entries. 
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Now we may prove the Theorem 1.2. 

Proof of the Theorem 1.2. Let us divide the above sum into sums over 
vectors J € V m , m > 1. Then, setting A = Ho, B = Ho — H and applying 
Theorem 2.3. we get formula (1.3), which completes the proof. □ 



3. The commutator method in dimension 2 

Theorem 1.2 is not valid for the Laplacian if d > 2. A way to overcome 
this difficulty for d = 2 is to consider differences of resolvents (this was 
pointed out to me by Y.Kannai). 

Lemma 3.1. Let R\ = (A — A) -1 be the resolvent of the Laplacian on a 
2-manifold M. Consider the difference of the resolvents R\$\ = R\ — R2\- 
This operator has a continuous kernel R\ t 2\(%i,%2) which has the following 
asymptotic representation on the diagonal as A — > oo 

oo 

flA,2A(*,aO~ £&»(*)(- (3-2) 

n=0 

where 

b n (x) = (n- 1)! • a n (x), n > 1. (3.3) 

Proof. The difference of the resolvents R\ t 2\ is a self-adjoint smoothing op- 
erator from L 2 (M) into the Sobolev space H 4 (M), and dimM < 4, therefore 
it has a continuous kernel (see [AK]). 

Let us now prove the expansion (3.2). 

Set z = —A, and let Rez > 0. We have (formally): 



oo 

f e -t(A +z ) dt= 1 

J z + A 

o 

Therefore 

7 1 °? °? °° 

i~ Zt dt. 



#(— ^r) = - / te-^ A+z Ut ~ - ftY t' 1+n a. 
dz y z + A J J J ^ 

n=0 

Let us integrate both parts from z to 2z: 

2z 00 oo 

2 a n t n e- zt dt. 



z + A 2z + A 
For n = we get: 



2 



n=0 



2z oo 2z 



J ao J e zt dtdz = J —dz = ao log 2 = a Q 
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For n > 1 we get: 

2z oo 2z 

t n e- zt dtdz 



f dz 1 
T(n + l)a n J = (n - l)!a n ( 



z n (2z) r 



z 

This completes the proof of the lemma. □ 

Lemma 3.4. Let R\(x\, x%) be the kernel of the resolvent R\ = (A — A) -1 
Then for A — > oo the difference of the resolvents has the following asymptotic 
representation on the diagonal: 

oo 

#a,2a(z,x) -p(x)- 1 X m (F™ + \x,x)-F™+ 1 (x,x)), (3.5) 

m=0 

where 

m /™\ 

. k a m—k 




X m = ^(-l) k [ l A fc A™- fe . (3.6) 



Proof. Setting A = H and applying Theorem 1.2 and Lemma 3.1 completes 
the proof. Note that we can define Xq = I without loss of consistence of the 
recurrent formula (2.4.). □ 
Now we are ready to prove an explicit formula for the heat invariants in 
two dimensions. 

Proof of the Theorem 1.4. Let (u, v) be conformal local coordinates in 
the neighborhood of the point x = (0, 0) (we can assume that x is the origin 
without loss of generality) on the 2-manifold M. As it was mentioned in the 
introduction this is always possible. The metric in such coordinates has the 
form 

ds 2 = p(u, v)(du 2 + dv 2 ), (3.7) 

and hence the Laplacian is given by 

A = -p(u,v)- 1 (d 2 /du 2 + d 2 /dv 2 ). (3.8) 

Our aim is to find the coefficient a n {x). In order to do this we have to 
collect all terms in the sum ( |3.5| )) containing (— A) _n . We use the following 
well-known formula (see [AK]): 

du ai dv a ' 2 [X,X> { A> 1 ' 4vr2 J (Ai(6,6) + l) m+ Ko^ 

where A x (£i,£2) denotes the symbol of the operator A x , and (—A) 2 m is 
the analytic branch of the power which is positive on the negative axis. 

By formula ( [3.S| ) we get \a\ = 2(|J|+r — n) = 2(m — n). Since the integral 
( |3.9[ ) vanishes unless both ct\ and 02 are even we may put a± = 2m — 2n — 2p, 
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«2 = 2p, where < p < 2(m — n). In polar coordinates this integral 
transforms into the product of two integrals: 



/ 



/-2m— 2n— 2p j-2p 



(A;(6,6) + i) m+1 



2tt 



= \j COs2m ~ 2n ~ 2P ^ siQ2P( P di P J 



(t/ P0 + l) m +! 





= p^ n+1 B{m -n-p + l/2,p+ l/2)B(m - n + 1, n), 

where po = p(0, 0), and B-functions stand for the values of the first and the 
second integral in the product respectively (cf. [GR], formulas 3.621(5) and 
3.194(3); note that formula 4.368(3) which could be used directly for compu- 
tation of the integral in (3.9) is inaccurate). Expressing the S-functions in 
terms of T-functions we find that (3.9) may be represented as (— \)~ n I mnp , 
where: 



T (-Um-nJ_ n m-n+l £(m - U - p + l/2)T(p + 1/2)1» 



4vr 2ru r(m+l) 

Since \a\ should be positive, we have m > n. On the other hand, as it was 
proved in ([AK]) the order of every operator [B,A;J\ is not greater than 
| J | +2r and hence |a| < | J| +2r which implies \J\ < 2n. Another statement 
from ([AK]) gives the same estimate on the length of the vector: r < In. 
Therefore if we are interested in the coefficient a n it is sufficient to consider 
X m with m < An. 

Together with (3.3), (3.5) and (3.6) this implies: 

m-n „,2m-2n-2p„,2p 

° ( p ^ mnP (2m -In- 2p)!(2p)! )|u = J° 10) 

Computing explicitly A™~ fc using (3.8) we may apply the obtained differ- 
ential operator to the last sum in (3.10). Afterwards we change for conve- 
niency one of the summation indices and omit the vanishing terms. Finally 
we obtain: 

4n m k—n 

*n{x)= E E Y, C nksmP(u,v) k - n A k (u 2k - 2n - 2s V 2s )\ u=0 , v=0 , 

m=n+l k=n+l s=0 (3-H) 

where the constants C n ^ sm are given by 

_ (-l) n ^ T(k + I - n - s + l/2)r(a + m-k -1 + 1/2) 
nksm ~ ^ k\l\{m- k-l)\{2k-2n-2s)\{2s)\ ^ ^ 
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This completes the proof of the Theorem 1.4. 



□ 



4. Curvature coordinates and an invariant expression for a n (x) 

IN A "GENERIC" CASE 

The commutator method essentially involves non-invariant coordinate 
representation. Nevertheless, using a special substitution it is possible to 
get the result in an invariant form, though in such a way that it is valid for 
in some sense "generic" germs of curvature functions on 2-manifolds. 

Definition 4.1. Let x be a point on the manifold M, and K(x) = Kq be 
the Gaussian curvature at this point. Let (z, w) be coordinates near the 
point x and suppose that in a neighborhood of x the map 

(z, w) - (z K ,w K ) = (K — K , AK - AK ), 

where AKq = (AK)(x), has a nonvanishing Jacobian. The coordinates 
{zk,wk) are called the curvature coordinates in the neighborhood of the 
point x. 

Lemma 4.2. The operator A$ = A' x may be written in curvature coordi- 
nates as 

A = -{Ed 2 /dz 2 K + 2Fd 2 /dz K dw K + Gd 2 /dw 2 K ), 

where 

2E = 2g u = -A(z 2 K ) = 2KAK - A(K 2 ); 

2F = 2g 12 = -A{z K w K ) = KA 2 K + (AK) 2 - A(KAK); 

2G = 2g 22 = -A{w 2 K ) = 2AKA 2 K - A(AK) 2 , 
and all the values are taken at the point x = (0, 0). 

Lemma 4.3. In the new coordinates 

u = V EG — F 2 zk v = Ewk — Fzk, 
the Riemannian metric {gij} takes the following form at the point x = (0,0): 

ds2 = E{ EG-F*) {du2 + dV ^ (44) 
where E, F , G are as in Lemma 4.2. 

Lemmas 4.2. and 4.3. are proved by easy computations. 

The proof of the Theorem 1.4. uses the conformal form of the metric in 
fact only at the point x = (0,0), therefore we may apply this theorem to 
(4.4), and expressing the result in curvature coordinates we obtain: 
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Theorem 4.5. Suppose that curvature coordinates exist in some neighbor- 
hood of the point x = (0,0) G M. Then the heat invariants a n (x) are given 
by the following formula 



where C n k sm are the same as in formula (3.12) 

Note that all terms in this expression are invariant — indeed, it is con- 
structed from the powers of the Laplacian and the Gaussian curvature taken 
at the given point. 

It is not always possible to introduce curvature coordinates — for ex- 
ample, on a standard round sphere the Gaussian curvature is constant 
and therefore such coordinates are degenerate. Nevertheless, they do ex- 
ist " generically" : 

Theorem 4.7. In the space J r of r- jets (r > 3) of germs of analytic func- 
tions K(u, v) at the point x = (0, 0), the set of r-jets for which the "curvature 
coordinates" are degenerate forms an algebraic set of codimension at least 
1. 

Proof. We demand the curvature function K(x), x € M be analytic near 
the point i £ M in some local coordinates (and hence in any). One may 
show that every germ of an analytic function maybe taken as a germ of some 
curvature function. Indeed, if we write the metric locally in conformal form 
(3.7) the Gaussian curvature is given by 



Clearly, for every analytic function K(u, v) there exists a function p(u, v) 
satisfying the above PDE with initial conditions K(0, 0) = Kq (by Cauchy- 
Kovalevskaya theorem — see, for example, [ES]). 

Consider the space J r of r-jets (r > 3) of germs of analytic functions 
K(u,v) at the point x = (0,0). Non-degeneracy of the "curvature coordi- 
nates" is equivalent to the nonvanishing condition for the determinant of the 
Jacobi matrix Jac(-fT, AK) at the point x. The equation Jac(-fT, AK) = is 
an algebraic equation in the space J r , therefore (since the equality is obvi- 
ously not identical) the set of its solutions is an algebraic set of codimension 
at least 1 (see [Br]) in the space J r . □ 



As a further check of the formula (3.11) we programmed it (together 
with its slight modification for diagonalized metrics) using the Mathematica 
software ([Wo]). 





(AK-AK y)\ x=m , (4.6) 




(4.8) 



Appendix 
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Here is the result of the work of the program for a\ (x) : 

( N _ 1 pj + Pi - PUUP ~ PwP 

ai[X) ~ 24vr p* 

On the other hand, using the formula (4.8) we get that 

„ pl + pl- puuP - PwP 

K = p ' 

and hence 

ai{x)= 1^ = 247' 
where r = 2K is the scalar curvature at the point x. This coincides with 
the well-known formula for a\{x) (cf. [G3]). 

Let us emphasize that our aim was just to test the formula (3.11); for 
computational purposes Mathematica is not an efficient tool for our method 

— it does not manage to calculate even 02 (x) for a general metric in reason- 
able time (though in particular cases of a "simple" conformal metric with 
p(u,v) = 1 / (ao + a\u + a2v) , or a spherical metric ds 2 = R 2 (d9 2 + sin 2 6dip 2 ) 

— here we used a modification of (3.11) for diagonalized metrics, — we found 
a,2(x), and the result also agreed with the already known (see [G1],[G3])). 
The main difficulty is to compute high powers of the Laplacian, which is a 
differential operator with non-constant coefficients, and therefore the num- 
ber of terms in the expression for A k grows exponentially in k. However, 
we think that programming the commutator method in a language like C or 
in some lower level language one may obtain completely explicit expressions 
for higher heat invariants and, hopefully, even for some unknown ones. 
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